Horizon in Random Matrix Theory, Hawking Radiation and Flow of Cold Atoms 
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We propose a Gaussian scalar field theory in a curved 2D metric with an event horizon as the 
low-energy effective theory for a weakly confined, invariant Random Matrix ensemble (RME). The 
presence of an event horizon naturally generates a bath of Hawking radiation, which introduces 
a finite temperature in the model in a non-trivial way. A similar mapping with a gravitational 
analogue model has been constructed for a Bose-Einstein condensate (BEC) pushed to flow at a 
velocity higher than its speed of sound, with Hawking radiation as sound waves propagating over the 
cold atoms. Our work suggests a three-fold connection between a moving BEC system, black-hole 
physics and unconventional RMEs with possible experimental applications. 

PACS numbers: 71.10.Pm, 72.15.Rn, 73.23.-b, 04.62.+V, 04.70.Dy, 03.75.Kk, 67.85. Hj 

Keywords: Random Matrix Theory, Luttinger Liquid, Hawking Radiation, Event Horizon, Cold Atoms, 
Bose-Einstein Condensate, QFT in Curved Space-time, Bosonization 



One of the most spectacular results in general rela- 
tivity is the prediction by Hawking [l], Q of a thermal 
radiation whose existence is only due to the presence of 
an event horizon near the black hole singularity. Unfor- 
tunately, direct experimental evidence of Hawking radi- 
ation is almost impossible (see, for instance, 0]). After 
the beautiful proposal by Unruh 0], there has been a 
growing interest in searching for similar phenomena in 
the so called analogous gravitational models One of 
such models can be realized in a stream of cold atoms 
with an abrupt change of the stream velocity v(x) at 
x = 0. In such system, the event horizon arises due to 
the fact that the phonon velocity c, which is larger than 
stream velocity for x < 0, abruptly gets smaller than 
v for x > 0. Thus the entire region x > is analo- 
gous to the black hole interior from where no particle 
can escape. However, unlike the interior of black holes, 
the region x > is accessible to measurements, in which 
case an experimentalist plays the role of a super-observer. 
In this realization, Hawking radiation is characterized by 
an intensity (or phonon density nix)) which is strongly 
anti- correlated not only at small distances \x — x'\ ~ n~ l 
but also at the mirror-reflected condition |ic + a:' | ~ n~ x , 
being exponentially small elsewhere Amazingly, ex- 
actly the same behavior was found for eigenvalue density 
correlations (two-level correlation function TLCF) in a 
certain exactly solvable random matrix ensemble (RME) 
0,1. 

In this letter we propose a connection between three 
very different areas of physics, namely Random Ma- 
trices, Quantum Gravity and Bose-Einstein Conden- 
sates (BEC). The relation between the latter two has 
been already established in @, Q. Here we focus on 
a non-standard invariant ensemble of Random Matri- 
ces described by a probability distribution P(H) = 
P(UHU- 1 ) of the form 



In |8j, the asymptotic behavior of the TLCF for this 
weakly confined ensemble was evaluated and it was found 
that, beside the usual translational invariant part already 
found in 0, there was an anomalous, non-translational 
invariant component of the cluster function: 
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which corresponds to repulsion between the level at a 
point x in the (unfolded) energy space and its mirror 
image at x' . 

The main result of this letter is that the TLCF with 
this anomalous component can be exactly reproduced 
within a scalar field theory (Luttinger liquid of levels (loj ) 
in a curved space-time with an event horizon. The same 
kind of theory has already been shown in Q to be a good 
low-energy approximation for the BEC system pushed to 
move faster than its speed of sound, as it was confirmed 
by subsequent numerical simulations The equiva- 
lence of the TLCFs for these theories strongly suggests 
some relation between them and unlock the possibility 
of extending techniques from one field to another. For 
instance, using the Luttinger Liquid description, we pre- 
dict an oscillatory component in the TLCF of 6|, which 
should be experimentally accessible for a BEC in the 
Tonks-Girardeau regime. 

Due to the interdisciplinary nature of these subjects, 
we first briefly remind the reader of the different ingre- 
dients needed to establish these connections, so not to 
hinder the physical understanding behind the formulae. 

Hawking Radiation in a BEC: When a BEC fluid is 
pushed to move faster than its speed of sound, a sonic 
black hole with its event horizon is generated. This 
horizon separates a region where elementary excitations 
(phonons) can propagate along the fluid flow or against 
it and a region where the upstream propagation is sup- 
pressed, because the fluid velocity is so high that it drags 
any sound wave with it. Physically, it is simpler to mod- 
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ulate the local sound velocity c, while keeping constant 
the density po and velocity vq of the condensate Q . One 
can then design the set-up to have c\ < v < c r , where C\ tI 
are the sound velocities to the left /right of the horizon, 
where c = vq. 

The fact that the boundary between the sub- and 
super-sonic regions of a fluid constitutes a horizon is rel- 
atively intuitive. Much less obvious 3 is the fact that 
the mere existence of a horizon kinematically implies the 
generation of Hawking radiation emanating from it 
Analogue gravitational models like this one are a topic of 
growing interest 0, since "table-top" experiments can 
be used to test predictions of quantum gravity and help 



in resolving some of its issues [111 ]. However, the experi- 



mental detection of Hawking radiation has been so far im- 
peded by its small temperature and intensity. The novel 
idea put forward in [6|] has been to investigate the TLCF 
as a signature of Hawking radiation, as pairs of phonons 
generated close to the horizon will remain highly entan- 
gled as they travel in opposite directions in the sub- and 
super-sonic region. 

Sound waves propagate according to a d'Alembertian 
(a Laplace-Beltrami equation in a curved metric). The 
fluid is 3-dimcnsional, but since the relevant degree of 
freedom is only the one along the fluid motion, two spa- 
tial dimensions can be traced out and the effective dy- 
namics is 1 + 1. In a flat space, the fundamental Green's 
function is G(u) cx — | In (u + u~), where u ± = t± J 
are the light-cone coordinates. The phase field for the 
condensate propagates according to this correlator, the 
number density field being related to it by a differential 
operator (see [6] for the details). 

The key element in [6] is that, when a sonic black-hole 
horizon is created, the outgoing modes are exponentially 
redshifted in frequency and they propagate according to 
the effective light-cone coordinate 
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(3) 



where k is a parameter known as the surface gravity at 
the horizon. The ingoing modes propagate normally ac- 
cording to u + . Taking Q into account for the propa- 
gator, the two-point density-density correlation function 
can be calculated, yielding a translational invariant as 
well as an anomalous contribution for x x' < 
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This correlator shows a characteristic off-diagonal peak 
connecting the sub- and super-sonic region of the fluid. 
In this low-energy, field-theory based prediction has 
been tested against microscopic, ab-initio numerical sim- 
ulation of the system showing an almost perfect agree- 
ment. Please notice that selecting c T> i — vq ± vq/2, (|U) 
correctly reproduces the denominator in 



The weakly confined RME: One exactly solvable RME 
with a confinement behaving like |T]) is 



V{E) = In [l + 2g n+1 cosh(2x) + q 
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(5) 



n=0 



with E = sinhx, q = e~ K and k > [|. The model has 
been solved using appropriate orthogonal polynomials (q- 
deformed Hcrmitc polynomials) and other methods and 
various quantities have been calculated in the literature. 
One can introduce a change of variable that makes the 
mean level density p(E) = Tr {5 (E — H)} constant, i.e. 
d£ 

(p(x)) = (p(E x )) —r^- = 1 . One of the most important 

Ua; 

peculiarities of this model is that the variable perform- 
ing this unfolding in the bulk of the spectrum behaves 
exponentially Q: 
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where A is an (uninfluential) constant of integration. 
The TLCF is defined as 



Y 2 (x,x r ) =5(x-x') 



(p(E x )p(E x ,)) - (p(E x ))(p(E x ,)) 
(p(E x )){p(E xl )) 



For p = e~ 27r l K <C 1 a semiclassical analysis can be 
carried out Q and one finds that the asymptotic behav- 
ior of the two-point correlation function is Y%{x, x') — 
Y£{x - x') 9{x x') + Y£{x, x') 9{- x x') , where 6{x) is 
the traditional step function. The normal translational 
invariant part of the correlation function is p} 



Y' 



X J = —r — 



sin 2 [7r(j: — x')\ 



4tt 2 sinh 2 [K(a; - x')/2] 



for x x > 0, (7) 



while the anomalous (non-translational invariant) part is 
given by ((2|) . As the anomalous component is the most in- 
teresting one, we do not concentrate on the normal com- 
ponent (O, even if it is also common to all theories con- 
sidered in this paper. The main difference between @ 
and (|1]) is the oscillatory part at the numerator of ^j, 
the derivation of which is the focus of the next section. 

RME and Luttinger liquid: If we interpret the en- 
ergy levels of a random matrix as the coordinates of a 
system of particles, the distribution of eigenvalues can 
be thought of as the equilibrium configuration of such 
a system of quantum interacting particles in ID, the 
phenomenon of level repulsion being encoded in their 
fcrmionic statistics and in the form of interaction. In 
this description, the probability distribution ([T]) becomes 
a confining potential for the particles. In condensed mat- 
ter, it is well known that the low-energy behavior of a ID 
system is not given by a Fermi Liquid theory, but instead 
by the so called Luttinger Liquid (LL). 
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The Luttinger Liquid paradigm stems from the fact 
that in ID a system responds as a whole to any per- 
turbation because particles cannot go around each oth- 
ers. Therefore, low-energy excitations have an intrinsic 
sound-wave nature and the system can be described in 
generality using only two parameters K and c within a 
quadratic action for a bosonic displacement field <&(x,i): 



5[#] = 
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- (9 T $) 2 
c 



c(d x $) 



(8) 

Contrary to Fermi Liquid theory, the Luttinger Liquid 
description is valid for systems with fermionic statistics 
as well as for bosonic systems, although the details of the 
mapping are of course different [l2| • In (J8J), c is the wave 
(sound) velocity (we will set h = c = 1 henceforth) and 
the parameter K encodes all the effective interactions of 
the original model: K — > oo for free bosons and decreases 
as the repulsive interaction increases, while for fermions 
K > 1 corresponds to attraction, K < 1 to repulsion and 
K = 1 means free fermions. In this low energy approxi- 
mation, the bosonic theory is quadratic and therefore the 
fundamental Green's functions is G(x) oc — | In (x 2 )- 

The particle density operator can be written in terms 
of the displacement operator as [l2| 

1 Ak 

p(x, t) = po - - d x <f> + cos [2tt Po x - 2$] + . . . , (9) 

where po is the background density - po = 1 in our RME 
after unfolding ([6]) - and higher harmonics terms arise 
from interference effect (Umklapp processes mixing left 
and right movers in the fermionic language). Ax is a 
constant that depends on the short-distance behavior. 

We propose that the low-energy effective theory for the 
invariant weakly confined RME defined by JT|) is given 
by the ground state of a LL in a curved metric, i.e. by a 
bosonic field whose dynamic is dictated by the action 



i 



2-kK 
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where g = det g^ v is the metric, i.e. ds 2 = g^d^d^ v . In 
this case the parameter k enters in the space-time metric 
g rather than directly as a temperature T like in ([5]) . As 
we are interested in the low-energy asymptotics of the 
theory, we cannot specify the metric uniquely, but its 
causal structure will suffice. 

In the construction, we are guided by the work done 
for the supersonic BEC [6]. We need a horizon, since 
this causes an exponential redshift in the modes, which 
leads to the appearance of Hawking radiation and of an 
anomalous two-point correlation function like ([2]). How- 
ever, the metric used in Q is physically motivated by the 
BEC setup, where time reversal invariance is explicitly 
broken, hence only the out-going modes are exponentially 
red-shifted. For the RME we will preserve T-reversal and 



therefore we will not have a radiation, but a bath with a 
finite temperature T. 

In [l(| , it has been shown that the traditional Luttinger 
liquid phenomenology (|8I9| with T oc k correctly repro- 
duces the TLCF ([7]) of a non-invariant RME, which is a 
one-parameter deformation of the classic Wigner-Dyson 
RME, leading to the translational- invariant result (0). It 
is remarkable that is common to both systems, while 
to obtain the anomalous part ((2|) we need to introduce 
the temperature in a non-trivial way through the metric 
(ITU]) and a horizon. 

Let us then start with the Rindler line element 



-y 2 dt 2 



— dy 2 = -y 2 du + du , 



(11) 



with u ± = t ± ^ln|y|, which is arguably the simplest 
1 + 1 metric with a horizon and it can be used as a 
model to capture the universal features of many phys- 
ical systems (for instance, as the dimensional reduc- 
tion of the Schwarzschild black-hole solution -in Kruskal 
coordinates- or as the effective metric observed by a uni- 
formly accelerated observer in flat background [13j). 

We can remove the singularity of the metric (fTTj) at 
the horizon y = by introducing a new set of null (light- 
cone) coordinates: 



±ie 



sgn(y), 



in terms of which we have a flat Minkowski space 



ds 2 



-du + du 



-dt 2 + dx 2 . 



(12) 



(13) 



where i2 ± = i ± x (the plus/minus refers to left /right 
moving modes, i.e. analytic/antianalytic sectors). 
Motivated by (|5I6[) . in (fTTj) we chose y = sinh(ra): 



ds 2 = - sinh 2 (ra) dt 2 + cosh 2 (ra) dx 2 



(14) 



x — — sinh(ra) cosh(Ki), t — — sinh(reir) sinhf/tt). (15) 

K K 

Lines of constant t are straight lines through the origin, 
while curves of constant x correspond to hyperbolae in 
the flat Minkowski space (|13|) . It is known that hyper- 
bolae are world lines of uniformly accelerated observers 
with acceleration a = a sinh _1 (Kcc) and indeed the set 
of coordinates (fT5|) is a modified version of those of the 
"traditional" accelerated observer [l]| and coincide with 
them far from the origin: 



e «|x| 

■r ~» — ; — cosh(/ct) sgn(x), t 



2k 



~2k~ 



■ sinh(Att) sgn(x) 



(16) 

In fact, the vacuum state of the Minkowski space is seen 
as a thermal state in (|14ll : an inertial observer would feel 
a black body (Hawking) radiation with temperature T = 
K/(27rfcs). This can be seen by noticing the periodicity 
in the imaginary time of (|16p . But it can be shown more 
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accurately by carefully decomposing the fields in modes 
in the two metrics and calculating their overlap [l3j |. 

We can now evaluate the density-density correlation 
function for the weakly confined RME using ©: 



Y 2 = 



1 



(d x ${x)d x ,${x')) 



(17) 



A 2 

-4cos(27r(x-a;'))(e 
2tt 



i2*(x) -i2*(x') 
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The non-oscillating part is simply 

K 



d x d x ,G(x,x') , (18) 



and for a Gaussian theory we have 

( e ia$(x) e -ia$(x')^ _ g - 3p { [*(x) -*(x')] ^) — e ~ T-^G{ X ,x') 
^ . (19) 

The (equal-time) Green's function in the coordinates 
given by (fL2")) - or (fTo) since we are interested in the 
asymptotics far from the origin- is 



G(x, x 



f |x|,|x'|3>l 
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X x' > 

«'<0 



Combining this with (fl7l4T9|) , it is straightforward to see 
that (HE]) are reproduced with K = 1 (and Ak = 1 [iH]). 

Conclusions: We have shown that the simplest 2D 
metric with horizons yields the anomalous two-point 
function ([2]) characterizing the weakly confined RME and 
which was also found for a supersonic BEC ^ . This sug- 
gests that all these theories share the same topological 
nature. This anomalous correlator is just another aspect 
of the emission of Hawking radiation at the horizons. 

Eq. (fl~2"j) and d3j) are at the heart of the equivalence be- 
tween these models. However, in ([3]) only outgoing modes 
are exponentially redshifted, since the BEC system is not 
time-reversal invariant (the horizon is created at a finite 
moment in time). This means that the Hawking radia- 
tion is an actual flow of particles coming from the hori- 
zon. Instead, the metric we constructed for the RME is 
time-reversal invariant and therefore there is no radia- 
tion, but a thermal bath in equilibrium (Hartle-Hawking 
effect p:'> ) . This is required since both modes have to be 
redshifted in (TT3|) in order to give the correct prefactor 
in (|2|) compared to ([7]). 

Two final remarks: (i) The Luttinger Liquid descrip- 
tion for ID systems (both fermionic and bosonic) directly 
leads to the emergence of oscillatory terms like those 
found in (|2l7p . In our model, since K = 1, these oscilla- 
tions decay as fast as the leading term. These contribu- 
tions were not present in the BEC result ([4]) [6] , since one 
needs to employ a proper regularization scheme, like the 
one provided by the Luttinger Liquid paradigm, to cal- 
culate next-to-leading corrections. However for K > 1 
they would be suppressed according to (flU)) . but for a 



cold atomic system in the Tonks-Girardeau regime these 
oscillations should be observable, (ii) Both the invariant 
RME ([I]) we studied here and the non-invariant, Power- 
Banded RME considered in seem to be described by 
a thermal Luttinger Liquid theory and share the transla- 
tional invariant TLCF ([7]). However, temperature in the 
invariant case is introduced in a non-trivial way through 
the Hawking effect, which leads to the anomalous com- 
ponent ((2]) as an indication of a non-trivial topology. It 
has been argued that the weakly confined invariant en- 
semble ([T]) spontaneously breaks its rotational invariance 
[1, IT3 an d would share many similarities with the system 
of [l0(. If this was the case, one could use the plethora 
of analytical techniques available for invariant models to 
study a non-invariant system. Our work provides new ev- 
idences of the connection between these two RMEs, but 
leaves many questions still unanswered. 

A microscopical derivation of the low-energy effective 
theory for the weakly confined RME could help in clari- 
fying this latter point and explaining the role of topology. 
Moreover, what the Hawking radiation represents should 
be understood in the original RME model as well as the 
interpretation of the time coordinate (and the tempera- 
ture) of the bosonic model, which in [1 01 ] was identified 
with the magnetic field threading the system. We should 
also remember that the RME |(SJ) is an exactly solvable 
system and therefore its interpretation as a gravitational 
analogue model has a lot to offer, for instance in address- 
ing the so called transplanckian problem. All these points 
will be the focus of our next efforts. 
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